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We develop a model-independent "box" parameterization of neutrino oscillations. Oscillation 
probabilities are linear in these new parameters, so measurements can straightforwardly determine 
the box parameters which can then be manipulated to yield magnitudes of mixing matrix elements. 
We examine the effects of unitarity on the box parameters and reduce the number of parameters 
to the minimum set. Using the box algebra, we show that CP- violation may be inferred from mea- 
surements of neutrino flavor mixing even when the oscillatory factor has averaged. The framework 
presented here will facilitate general analyses of neutrino oscillations among n > 3 flavors. 



I. INTRODUCTION 

If neutrinos have mass, then the states of definite mass, i^, may be distinct from the states of definite flavor, v a . 
(We will use latin indices i, j, k, ... to refer to mass states, and greek indices a, f3, 7, . . . to refer to flavor states.) If 
the eigenvalues of the neutrino mass matrix are non-degenerate, then neutrinos may change flavors, or oscillate, as 
they propagate. The long-standing solar neutrino deficit [Q, the atmospheric neutrino anomaly [0, and the recent 
results from the LSND experiment || can all be understood in terms of oscillations between neutrinos. Traditional 
oscillation analyses generally assume no more than three generations of neutrinos. But resonant oscillations for the 
sun, oscillations for the atmosphere, and the LSND data each require a different neutrino mass-squared difference in 
neutrino oscillations are to account for all features of the data Q . Since three-neutrino models can have at most two 
independent mass-squared differences, a sterile neutrino is apparently needed to reconcile all the data while retaining 
consistency with LEP measurements of Z — > vv Several four-neutrino analyses appear in the literature [Q|3). 
If, however, statistical or systematic errors in the data evolve in the future, or if some data turns out to have an 
explanation other than neutrino oscillations, then three-neutrino oscillations may be sufficient. One recent analysis 
of the complete set of data in the three-neutrino framework is given in reference . 

Oscillation probabilities depend on products of four mixing-matrix elements. Several parameterizations of the mix- 



ing matrix in terms of rotation angles have been introduced, beginning with the pioneering work of Kobayashi and 
Maskawa B]. When three or more neutrino generations are included in the oscillations analysis, the oscillation prob- 
abilities become complicated functions of the neutrino mixing angles. But oscillations are observable and therefore 
parameterization-invariant. One must ask if there is not a better description of oscillations which avoids the arbi- 
trariness of angular-parameterization schemes. In this paper, we introduce a new "box" parameterization of neutrino 
mixing valid for any number of neutrino generations. Oscillation probabilities are linear in the boxes, enabling a 
straighforward analysis of oscillation data. In what follows, we develop the algebra of the boxes and the unitarity 
constraints on that algebra. To conclude we illustrate the boxes' reduction to a basis in the case of three generations. 
A phenomenological analysis of existing oscillation data will be performed in a future publication. 



II. THE STANDARD FORMULATION OF NEUTRINO OSCILLATIONS 



The probability for a neutrino to oscillate from v a to up is given by the square of the transition amplitude: 
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where n is the number of neutrino generations and 



For relativistic neutrinos, <$>ij is given by 
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With a little bit of algebra, the oscillation probability may be brought into the form 
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The probability for an antineutrino to oscillate from u a to Vp is obtained from P [x) by replacing V with T^*. This 
is equivalent to changing the sign of <f>y . P (x) therefore differs from P (x) merely in the sign of the second term 

"alpha -» VjS "a, -> ^ 

in equation (^) or of the exponent in equation (0). Note that P (ac) = P (a;) , as required by CPT-invariance. 

v a -> ^ -» v alpha 



Let us first make contact with the familiar case of two neutrino flavors, n = 2. An arbitrary 2x2 unitary matrix 
will have one rotation angle parameterizing the real degree of freedom, and three phases. But all three phases may be 
absorbed into the definitions of Dirac fermion fields (or will cancel in oscillation probabilities for Majorana neutrinos 
0,0), so a 2 x 2 mixing matrix V has the simple form of a rotation matrix: 
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Since the matrix is explicitly real, the term in equation (|5|) involving imaginary matrix elements disappears, and the 
oscillation probability in the two-flavor case is 
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This result is simple, and the mixing-angle parameterization is a natural choice in the two-flavor situation. 

The math becomes more complicated when we include three generations. An arbitrary 3x3 unitary matrix has 
three real degrees of freedom and six phases. 2n — 1 = 5 phases are not observable [ pTj|JT^ ], since this is the number of 
phase differences among the fermion fields which may be absorbed into field redefinitions. The three-generation mixing 
matrix is therefore described by three mixing angles and one phase. The original choice of these four parameters, by 
Kobayashi and Maskawa to describe quark mixing, is perhaps the best known parameterization. Their choice of V is 
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where c a = cos# a , and s a = sin6> a . We will refer to this particular choice of angles as the "standard" or "KM" 
parameterization. 

In the standard KM parameterization, the phase only appears in the lower right-hand sub-block of the matrix. 
There is arbitrariness associated with the placement of the phase, since we absorb five relative phases into the field 
definitions. Clearly the location of the phase cannot be measurable. Indeed, the CP-violating effects of the phase are 
contained in a single function of the phase called the Jarlskog invariant J, which in the standard KM parameterization 
has the form |13)| J — c\s\c 2 s 2 c 3 s 3 sin 5. Because of the arbitrariness of the phase convention, the phases of individual 
matrix elements are not observable; only the magnitudes of matrix elements are observable. 
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The observable oscillation probabilities are quite complicated functions of the mixing angles in angle-based param- 
eterizations. As an example, consider the product ^22^2*3 V32V33 appearing in the — > v T oscillation probability. In 
the standard KM mixing matrix, this product is given by 

^22^*3^2^33 = 44 [44(4 + $4 + ?4 COS 2(5) - C?] 

+^ (1 + 4) (4 -sl)(s 2 3 -4) cot S + n = 3. (9) 
s i 

Besides being ugly, thereby motivating the often-made two-generation approximation, the expression (but not its 
value) on the right-hand side of equation (^) is convention-dependent. Other parameterizations of the unitary mixing 
matrix involving three angles and a phase are equally valid and yield similarly complex oscillations probabilities. Our 
development of a model-independent parameterization has been motivated by the arbitrariness and complexity of the 
traditional approach. 



III. A NEW PARAMETERIZATION 
A. Boxes Defined, and Their Symmetries 

The immeasurability of the mixing matrix elements in the quark sector has been addressed by numerous authors, 
such as those of references ]l^~[r^|. Measurable quantities include only the magnitudes of mixing matrix elements, the 
products of four mixing-matrix elements appearing in the oscillation probabilities, and particular higher-order func- 
tions of mixing-matrix elements [ ]l5| , |l8| |. As evidenced in equations (Q), (J|), and (||), neutrino oscillation probabilities 
depend linearly on the fourth-order objects, 

ai n 0j ee V m V^V Pj vl = VaiV^-V^j, (10) 

which we call "boxes" since each contains as factors the corners of a submatrix, or "box," of the mixing matrix. For 
example, the upper left 2x2 submatrix elements produce the box 

n o 22 = v 11 v 1 * 2 v 2 \v 22 . (11) 

These boxes are the neutrino equivalent of the "plaques" used by Bjorken and Dunietz for another purpose in the quark 
sector in reference fl6| . The name "box" also seems appropriate in light of the Feynman diagram which describes the 
oscillation process (Figure |l|) . 
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FIG. 1. Interference diagram contributing to the box which appears in the — > oscillation probability, as the 

neutrino travels a distance D. Arrows on the diagram represent fermion number, not momentum. 

In general, each of the box indices i, j, a, and (3 may be any number between 1 and n, the number of neutrino flavors. 
We therefore initially have n 4 possible boxes. Examination of equation (|l0|), however, reveals a few symmetries in the 
indexing: 
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If the order of either set of indices is reversed {id est, j <-> i or (3 <-» a), the box turns into its complex conjugate; 
if both sets of indices are reversed, the box returns to its original value Km. And if V is replaced by V\ then 
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B. Degenerate and Nondegenerate Boxes 

Equation ([l2]) demonstrates that boxes with a = or i = j, are real. Indeed, these are given from equation ([h]) as 

(13) 

Those boxes with both sets of indices equal are 



l a aj = \V ai \ 2 \V aj \ 2 , and ai U m = |V Qi | 2 |%| 2 . 
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We call boxes with one and two repeated indices "singly-degenerate" and "doubly-degenerate," respectively. Boxes 
with q / /J and i ^ j are called nondegenerate. As can be seen from equation (|^), singly-degenerate boxes with 



repeated flavor indices enter into the formulae for flavor-conserving survival probabilities, but not for flavor-changing 
transition probabilities. Degenerate boxes with repeated mass indices (including the doubly-degenerate boxes) do not 
appear in any oscillation formula. Degenerate boxes may be expressed in terms of the nondegenerate boxes, as will 
be shown shortly. This possibility and the symmetries expressed in equation ( |l2] ) allow us to express combinations of 
boxes in terms of only the nondegenerate "ordered" boxes for which a < (3 and i < j. 

The number of flavor- index pairs satisfying a < (3 ordering is N = (^j = "^" 2 ^ . N mass-index pairs similarly 
satisfy i < j ordering. Thus, the number of ordered nondegenerate boxes is N 2 . The number of doubly-degenerate 
boxes is n 2 . The number of ordered singly-degenerate boxes is \v? [n — 1) = nN for flavor-degenerate boxes, and the 
same for mass-degenerate boxes, yielding 2nN total singly-degenerate ordered boxes. For n — 3, N — 3 too, and we 
have nine ordered nondegenerate boxes, nine ordered singly-degenerate boxes, and nine doubly-degenerate boxes. For 
n > 3, matrix elements V a i with a > 3 are not easily accessible since a > 3 labels heavy charged leptons beyond the 
usual e, (i, t. Therefore, in the absence of processes with sufficient energy to excite the heavy charged lepton degree of 
freedom, boxes with flavor indices exceeding 3 are not physically accessible either. The number of accessible ordered 
nondegenerate boxes is 37V. The number of accessible ordered singly-degenerate boxes is 37V with flavor degeneracy 
and 3n with mass degeneracy, yielding 3(n + N) for the total. These counts are recapped later in Table Q, along with 
other box counts. 



C. Oscillation Probabilities in Terms of Boxes 



Using the symmetries expressed in equation (12), equation (g) becomes 
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where we have defined the shorthand 




When individual oscillations 



cannot be resolved due to uncertainties in energy or position, the sin 2 ^ and sin 2$^ are effectively averaged, giving 



the transition probabilities 
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These probabilities are equivalently given by 
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which may be a more familiar expression; the equivalence of these expressions will be shown later in this paper. 
The survival probabilities P (x) may be found from the transition probabilities P (x) by 

u a — y u a u a — > Ur 

p (x) = i - 53 p (*) ■ ( 18 ) 

Survival probabilities are more simply expressed in terms of degenerate boxes, or |V|s, rather than nondegenerate 
boxes. They are 



P (x) =l-4^)£ Sin2 ^ = 1 ~ 4 Y, Hl^iH^I 2 sin 2 (19) 
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If oscillations are averaged, the survival probability is 
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The third equality here results from unitarity of V. 

Interchanging a <-> (3 in equation (JTq) gives the time-reversed reactions P (x) : 



P (x) = -2 V V [2 ^ w sin 2 $y + <f m sin 2$*-,] + <5 Q/3 , (21) 

p i=l j>i 

so a measure of T-violation, or equivalently CP-violation, in the neutrino sector is 
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P (x) -P (x) =P (x) -P (x) =4V V T w sin2$y. (22) 



Ignoring possible CP- violating phases in the mixing matrix, there are = N real parameters determining 

V. Determining these N parameters determines the complete mixing matrix. Conveniently, there are N transition 
probabilities P (x) — P (x) . Thus, all of the information in the mixing matrix is contained in the N transition 

V a — > Va Up — > v a 

probabilities. In this sense, they form a convenient basis for determining all oscillation parameters. Of course, if the 
same transition probability is measured at two or more different distances, then all TV transition probabilities may 
not be needed to determine V. 

Allowing CP-violation in the mixing matrix, there are N real parameters and h(n — l)(n — 2) phases, for a total 
of [n — l) 2 parameters. With CP-violation, however, there are 2N — n(n — 1) independent transition probabilities 
P (x) . The number of transition probabilities exceeds the number of independent parameters, so they again form a 
convenient basis for determining the mixing matrix. In reality, only the three flavor indices e, fx, r are easily accessible, 
as alluded to at the end of the previous section. Moreover, some of the N parameters in the mixing matrix, namely 
those which rotate sterile states for n > 5, are not accessible. The inaccessibility issue complicates the counting. We 
do not pursue it further here. 
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The transition probabilities for which a ^ (3 in equation (15) may be expressed in matrix form. The matrix of 
boxes will necessarily be an N x N matrix since the vector on the left-hand side and the sum on the right-hand side 
of equation ( |l5| ) each contain N components or terms. For three flavors, we have 

\ 



V{n = 3) = 





(x) 


P 


(x) 

— > v T 


u 


(x) 
— > y T 



= -4 Re(B) 5 2 ($) + 2 Im(B) S(2$), 



(23) 



where 



B = 





e2 ^ 3 


el n M 3 


^□r2 




^□r3 


el Dr2 


e2 n 


Cl Dr3 



^ sin fe $12 



and S*($) = 



sin fc $ 



2:i 



sin $ 



13 



n = 3. 



(24) 



For the time-reversed channels, or for the antincutrino channels, the sign of the Im(£>) term is reversed. The definitions 
of V, B and S k ($>) may be extended to any number of flavors. Adding a fourth, perhaps sterile, neutrino flavor increases 
the dimensions of B to 6 x 6, since there are now six independent transition probabilities and six mass differences. The 
number of ordered nondegenerate boxes increases from nine to thirty-six, but only eighteen of them are accessible. 
For the six flavors of neutrinos in mirror-symmetric schemes, B is a fifteen- by-fifteen matrix, and we have two hundred 
twenty-five ordered nondegenerate boxes. But the number of accessible boxes is just forty-five. 

Our parameterization is especially well-suited for considering higher numbers of generations. While using boxes in 
n > 3 situations may be cumbersome, it is much less so than extending a mixing-angle parameterization to higher 
generations. Our matrix B of boxes merely acquires extra columns when new flavors are introduced; extra rows are 
not accessible at energies below new charged-lepton thresholds. Furthermore, oscillation probabilities are linear in 
boxes, no matter how many generations. 



D. A Return to Two Generations 

The boxes may be used to illustrate when the two-flavor oscillation formula (H) is a meaningful approximation and 
when it is not. For this, we assume ordering of the mass indices mi < W2 < • ■ ■ < m n . If the transit distance L of the 
neutrino is just large enough to impact upon the shortest oscillations length , but not upon the next-shortest, 

then 

sm 2 2d e J/ = -4% n (25) 
8 



is meaningful. If one amplitude lipj among those associated with Amy satisfying Am^-^^J 2 dominates the others, 
then 

Bk 2 20^ = -4« w (26) 
is a meaningful approximation. If there is a dominant mass m„ 3> m n -i, ■ • ■ , mi, then for -ffi&L -C Am 2 P — , 

n-l 

sin 2 20^/ = -4 ^ (27) 
i=l 

is meaningful. We will show later (equation (pPf)) that the right-hand side of equation ( p7j ) may be written as 
4 |V QIl | 2 |V/3„| 2 . If the largest masses are nearly degenerate, with degeneracy q, then the above becomes 

n—q n 

sin 2 2^/ = -4£ £ <%, (28) 

z— 1 j—n — q+1 

again meaningful. However, if two or more well-separated Am^-s exceed then the two-flavor formula is inappli- 
cable. 



E. Relations between Mixing Matrix Elements and Boxes 

Neutrino oscillation experiments will directly measure the boxes in equation ( |l5| ) , not the individual mixing matrix 
elements, V a i. But one would like to obtain the fundamental V a i from the measured boxes. We develop here an 
algebra relating boxes and mixing matrix elements. 

Some tautologous relationships between the degenerate and nondegenerate boxes are easily confirmed using equation 
(|l0|); they hold for any number of generations: 

\V ai f\V aj \ 2 = ai n aj = ain %" tnXj , A + a), (29) 

= ai u 0i = aia tJ n0V , (x*y*i), and (30) 
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In these equations and what follows, a, /?, 7, i, j, and k will usually be reserved for indices that are chosen at the 
start of a calculation, while other indices such as x, y, rj, and A primarily represent "dummy" indices which are chosen 
arbitrarily except to respect the inequality constraints following equations such as equation (^l|). 



The tautologies ( p9| ) to (|3l]) become evident not only by algebraically using the definitions of the boxes, but also by 
considering a graphical representation we have developed. This representation is discussed in detail in Appendix |A|, 
and it proves quite useful when considering relationships between boxes. Figure || illustrates the representations of 
"□22 and ( n n 2 2 
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FIG. 2. The graphical representation for (a) the box CI22 , and (b) its inverse ( O22 J ■ Vertical arrows point from the 
matrix elements which are not complex conjugated in the box to the complex-conjugated elements. While shown here in three 
generations, this representation may be extended to higher generations merely by adding more dots. 

To illustrate an example of equation (ps|), the degenerate box 12 Di3 is "created" in Figure ||a when one uncanceled 
vertical arrow enters and another leaves each of the matrix elements V12 and V13, while arrows at all other points 
cancel. Those arrows are the result of the combination of ordered boxes given in (|2^), as shown in Figure [|b. 
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FIG. 3. The graphical representation in three generations for (a) the degenerate box 12 Qi3 , and (b) the combination of 
nondegenerate boxes which "create" it. 



Equations (|29|) and (|30| ) are themselves special cases of the more general 
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and the analogous relation 
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(33) 



These relations hold for both degenerate boxes and nondegenerate boxes. Including disordered boxes adds no new 
information, so we will as usual consider only ordered boxes. 



This index rearrangement of equations (|32|) and (33) may straightforwardly be extended to obtain higher-order 
relationships (trilinear, quadrilinear, et cetera) in the boxes. These higher-order relationships do not enter into 
oscillation probabilities, so we will not examine them here. 

We may express \V a i\ = ( ai ^ai ) 4 in terms of three singly-degenerate boxes by setting a — (3 in equation (|3C 
Then, using equation (|29| ) to substitute for the singly-degenerate boxes yields an expression for the doubly-degenerate 
box in terms of nine nondegenerate boxes: 



on □ an □ 



ax u Ti ™u ax <*vo pi ™n (y "*d w 
f*n Cv a vn ux ax o^ y 



t =/= a =/= a 

C ¥= P + <* 
/i 7^ w 7^ a 



(34) 



As noted earlier, flavor-degenerate boxes give the probabilities for flavor-conserving oscillations measured in disap- 
pearance experiments, while nondegenerate boxes give the probabilities for flavor-changing oscillations measured in 
appearance experiments. Since the degenerate and nondegenerate boxes are related, individual \V a %\ may be deduced 
from an appropriate set of measurements of cither kind. 

One can obtain a relationship similar to equation ( |34"| ) by setting i = j in equation ( p9[ ) and then using equation 
( ^p| ) to substitute for the singly-degenerate boxes: 
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In the three-generation case, equations (Q) and (|3^) are uniquely specified by the index constraints in brackets, 
and they are equivalent to each other since all nine boxes are used in both equations. For example, 
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with a = i = 1 holds with any number of generations, but it is the unique 5 on 4 box representation of |Vn| 4 in three 
generations. Another example in three generations is 



V21I 4 = 



21 n Hn 21 n 1:L n 12 n 

u 32 u 22 u 33 u 23 u 33 



HD 32 "033 22 D33 12 D23 ' 

again true for all n and unique in n = 3. The equalities in equations (pi) and ([37]) are illustrated in Figure 0. 
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(a) (b) 

FIG. 4. The graphical representation for (a) | VT 1 1 4 , and (b) |V2i| 4 . Notice that each matrix point except the one being 
represented has an equal number of vertical arrows as horizontal arrows entering and leaving. The element being represented 
has two vertical arrows leaving (representing V^) and two entering (representing V*?) to produce |F Q i| 2 . Appendix |a] provides 
more details. 

If n > 3, equations (|34"| ) and ( |35|) still hold, but the indices of each are not specified uniquely. In higher generations 
the indices may take on values larger than 3; if n = 6, equation (Q) may contain up to six different flavor indices but 
only three mass indices, while equation ( |35| ) may contain more mass indices and only three flavor indices. 

Before continuing, we should note that all of the relationships in this section follow from the definitions of the boxes 
in equation (0) and so are valid for any matrix, unitary or otherwise. The constraints of unitarity, to which we will 
turn our attention after we count the number of independent boxes, will provide us with expressions for |V^i| 4 (found 
in equations ( |68| ) and (p9|)) which are easier to manage than the expressions (|34|) and ( |35| ) above. 

F. The Number of Independent Boxes 

An n x n arbitrary matrix has n 2 elements, which may be complex, leading to 2n 2 parameters. The neutrino mixing 
matrix, however, is unitary, and has n 2 unitarity conditions as constraints. Of the n 2 remaining parameters, ^n(n— 1) 
(the number of parameters in a real unitary, or orthogonal, matrix) may be taken to be real rotation angles, and the 
remaining ^n(n + 1) are phases making the matrix complex. We may redefine the 2n neutrino and charged lepton 
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fields to remove 2n — 1 relative phases from the mixing matrix, leaving |n(n + 1) — (2n — 1) = i(n— l)(n — 2) phases 
measurable in oscillations. The number of independent real and imaginary box parameters should be the same as 
the number of real and imaginary CKM parameters [ fl8| |. As we will see, unitarity relates the imaginary and real 
box parameters, so one may choose more real boxes and fewer imaginary boxes as the independent basis of (n — l) 2 
elements. In fact, one may choose a basis of only real boxes. For n > 3, not all of the boxes are generally accessible 
(and not all of the (n — l) 2 mixing matrix parameters are accessible), so the basis of boxes may be smaller. The 
parameter counts for the mixing matrix and the boxes are summarized in Table Q. 
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TABLE I. Parameter counting for the mixing matrix and boxes. N is shorthand for the frequently-appearing combination 



[j = ^n(n — 1). n/ represents the number of active neutrino flavors, not to exceed 3, and Nf = 



Number of generations 


n 


2 


3 


4 


6 


Params. for arbitrary matrix 


2n 2 


8 


18 


32 


72 


Unitarity constraints 


n 2 


4 


9 


16 


36 


Relative phases 


2n - 1 


3 


5 


7 


11 


Real params for unitary V 


N 


1 


3 


6 


15 


Remaining phases in V 


I(„-l)(n-2) 





1 


3 


10 


Physical params in V 


(n - l) 2 


1 


4 


9 


25 


Initial boxes 


n 4 


16 


81 


256 


1296 


Doubly-degenerate boxes 


n 2 


4 


9 


16 


36 


Ordered singly-degenerate boxes 


2nN 


4 


18 


48 


180 


Ordered nondegenerate boxes 


N 2 


1 


9 


36 


225 


Independent Re( c "D ( 3j )s 


N 


1 


3 


6 


15 


Independent Im( al ^pj )s 


i(n-l)(n-2) 





1 


3 


10 


Accessible boxes 


2 2 

n f n 


16 


81 


144 


324 


Accessible doubly-deg. boxes 


n f n 


4 


9 


12 


18 


Accessible ordered singly-deg. boxes 


rifN + Nfn 


4 


18 


30 


63 


Accessible ordered nondeg. boxes 


N f N 


1 


9 


18 


45 
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IV. UNITARITY RELATIONS AMONG THE BOXES 



Unitarity requires that 

n 

Y^VniVjj = By, and (38) 

•7=1 
n 

Y^V ay V; y = 5 aP . (39) 

These equations are not independent sets of constraints at this point: equation (|3^) states V^V =4, and equation ( |39| ) 
states =4. We can, however, use these equivalent equations to obtain two separate sets of constraints on boxes by 
multiplying equation ( |38"| ) by V^Vxj and equation ( p9| ) by V* x Vp x : 



53 = V Xi Uxi Sij, and (40) 

r,= l 



J2 ^Opx = V <W- (41) 



Isolating the manifestly degenerate boxes from the nondegenerate boxes, equation ( f40| ) becomes 



53 "*D Aj = V Ai DAi ~ , (42) 



and equation ( |4l| ) becomes 



53 "VQ^ = V M Q M <W - • (43) 

One can show that equations (^2j) and p^), although distinct, are related by the symmetry that takes al ^pj — > 
la ^*p , which we identified earlier as the symmetry — ► V*. We will find several pairs of equations related by this 
symmetry which differ only in whether the sum is over a flavor index or a mass index. 
Summing equation ( ppj ) over A in the i ^ j case, we find 

n n n n 

= E E = E Al °^ + 2 E E ( 44 ) 

A=l77=l A=l A=1?)<A 

Comparison of equation with equations ( p4|) and (|29|) reveals an interesting property of the matrix B: 

n 

53 Be(B) = -^\V Xi \ 2 \V Xj \ 2 , (45) 
columno/B A=l 

where the sum is over a column of B specified by fixed i and j. We may in a similar manner show that 

1 " 

53 Re(B) = --Y,\V ax \ 2 \V 0x \ 2 , (46) 
rowo/23 x=l 
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with constant a and (3 specifying the row of B to be summed. Sums over rows are preferred to sums over columns 
phenomenologically when n > 3, because complete rows are accessible whereas complete columns are not. 

For n > 3, some Rs will be positive while others are negative. An upper bound on the values of the positive Rs may 
be established from the Schwarz inequality for two-component complex "vectors" V a = (V a i, V a j) and Vp = (Vpi, Vpj): 



\v ai \ 2 + \v aj \ 2 



\V(K\ 2 + \V Pj \ 2 



(47) 



A bit of algebra gives the desired result: 



\v ai \ 2 \v^\ 2 + \v aj \ 2 \v 0i \ 2 



(48) 



Later we will actually derive an equality in three generations for "Hpj in terms of the four |V|s appearing in this 
equation (found in equation (p7|)). 

An alternative way to obtain constraints from unitarity is to start with the definition of the boxes (^) and use the 
unitarity of the mixing matrix: 



EW 

ri=jta 



E^ 

A#/3 



(49) 



After a bit of algebra, this becomes]] 



EE = -%(-i+i^i 2 + i^i 2 ) 

r,jta A^/3 



(51) 



The number of unitarity constraints from equations (40) and (|H|) , both real and imaginary, grows as n 3 , while the 
number of ordered boxes grows as n 4 , so additional relationships between ordered boxes must exist to identify the 
(n — l) 2 independent Js and Rs. These additional identities will come from the definitions of the boxes, and the 
resulting tautologies (|2^) and (|30"|), each of which grows in number as n 6 , implying a high degree of redundancy. Some 
of these relationships are explored below in equations ( |74] ) to 

1 This relation also follows from equation (p{i|): 



EE% =E E" <[ 



J Aj 



Ed^i 2 ^ 



J Aj 



A#/3^a 



\^f3 \n=l 



A#/3 
D /3j 



= ^Ei Va< i 2 - E ai °» - 

A^/3 \A=1 

= Sij (l — V^] 2 ) — \V a i\ 2 Sij + a! D/3j 



(50) 
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The constraints (|4(j) and ([H]) hold independently for the real and imaginary parts of each sum. We will first 
explore the implications of these constraints for the imaginary parts of boxes, before turning to the implications of 
these constraints for the real parts. 



A. Unitarity Constraints on the Imaginary Parts of Boxes 



Because the right-hand sides of the equations are manifestly real, as are terms on the left-hand side involving 
degenerate boxes, the sums of nondegenerate boxes in equations (^) and ([II]) must be real. This leads to imaginary 
constraints of the formpl 



J2 =0, and 

E 



= 0. 



(54) 
(55) 

In three generations, each sum in equation ( |54| ) or (^5|) contains only two terms, leading to the equality (up to a 
sign) of two Js. For example, choosing A = 1 in three generations yields 



%i = -%i, n = 3 



(56) 



so 



22 



32- 



n = 3. 



(57) 



We may continue in a similar manner and show that every J is related to J = ¥22 JL3|] , giving 



' J J -J 



Im(S) 



n = 3. 



(58) 



J J -J 
\-J -J J J 

One consequence of the equality up to a sign of all Js in three generations is that if any one V a i is zero, then all Vpj 
id there can be no CP- violation. 



The constraints (^) and (|55| ) may be written exclusively in terms of ordered boxes as 

7?<A T)>A 



(52) 
(53) 



While expressing everything in terms of ordered boxes is our goal, we will find it more convenient to use the complete sums of 
equations (^) and (JH^) in mathematical manipulations and switch to ordered boxes at the end rather than to deal with the 
two separate terms in equations (152) and (p3D. 
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For n > 3 generations, the sums in equations ( p4| ) and ( pq ) contain more than two terms, so individual Js become 
equal to sums of other Js rather than just equal to another individual J. In four generations, we know from counting 
phases in the mixing matrix, or more rigorously from the work in reference ]T^| that we should find three independent 
Js out of the 36 ordered Js. This reduction, unlike the reduction in three generations, is not achieved solely by 
equations of the forms ( |54| ) and (|55|). No more than 27 independent constraints may be derived from those equations, 
so at least 6 constraints on the imaginary parts of boxes in four generations will come from the expressions derived 
below which relate real and imaginary parameters. 



B. Unitarity Constraints on the Real Parts of Boxes 

We now consider the real parts of the constraints ( ^0| ) and (|il]) , 

considering first the homogeneous constraints for which the Kronecker delta is zero. These relationships give the 
singly-degenerate boxes as sums of ordered boxes: 

\V M \ 2 \V Xj \ 2 = A *D AJ = - £ *U Xj =-J2 *Xj, i * J, and (59) 

\V ax \ 2 \V 0x \ 2 = ax U 0x = - J2 ay V Px =-J2 °K/9«, a ^(3. (60) 

y^x y^x 

For three generations, each of the sums contains two terms, allowing us to express the singly-degenerate boxes in 
terms of two nondegenerate boxes measurable in neutrino appearance oscillation experiments. 

For fixed in equation (p39|), A can take n possible values, implying n constraint equations. N ordered non- 

degenerate boxes appear in these n equations. Thus, for N < n, which is true for n < 3, the unitarity constraint 
( p9| ) and (|6^) may be inverted to find a nondegenerate box in terms of singly-degenerate boxes. Adding the sums in 
equation (^30|) for x = i and x — j, then subtracting the x = k sum yields the desired expression: 

^i = -^(|K i | 2 |%| 2 + |K J -| 2 |%| 2 -|y afe | 2 |^| 2 ), n = 3, (61) 

where we have used the property that Rs are real and therefore unchanged under mass (or flavor) index interchanges. 
Many sources, such as reference ||, use these sums as the coefficients of the oscillatory terms in three-flavor oscillation 
probabilities to write 

P (x) =2(sin 2 $ 12 +sin 2 $ 1 3-sin 2 $ 23 )|14i| 2 |l/ /3 i| 2 

V a -> Vfj 

+ 2 (sin 2 $ 12 - sin 2 $ i3 + sin 2 $ 23 ) \V a2 \ 2 \V p2 \ 2 (62) 
+ 2 (- sin 2 $ 12 + sin 2 $ 13 + sin 2 $ 23 ) |V Q3 | 2 |Va 3 | 2 , n = 3. 
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Similar manipulation of equation ( |60| ) gives an alternate expression in term of the flavor triad (a, (3, 7): 

%i = ~ (\V ai \ 2 \V a3 \ 2 + \V 0l \ 2 \V M \ 2 - |V 7 i| 2 |V 7J f) , n = 3. 



(63) 



A simpler derivation of equations (pj) and (p3J) is available. Three-generation unitarity gives 



-v ak v, 



(64) 



Squaring both sides of this equation then leads directly to equation (pi]). Generalizing this derivation to n > 3, we 
see why a single nondegenerate box cannot be expressed as a sum of degenerate boxes for n > 3, but we also obtain 
some new relations among boxes. For example, with n — 4, we may square both sides of 



v al vz + v a2 vL = -VaiVZ - v al v; x 



(65) 



to obtain 



°ft/32 — °K/34 = 



1 



|K 3 | 2 \v^\ 2 + \v ai \ 2 \v Pi \ 2 - \v a i\ 2 |^i| 2 - \v a2 \ 2 \vp 2 \ 2 



n = 4. 



(66) 



It has been shown that knowledge of four |V|s completely specifies the three-generation mixing matrix, provided 
no more than two |V|s are taken from the same row or same column |l9[] . Therefore, it may prove useful to use 
three-generation unitarity to re-write equation (p3[) in terms of just four |V|s. The result is 



1 



1 - \v m \ 2 - \V a3 \ 2 - |^| 2 - \V PJ \ 2 + \V m \ 2 |%| 2 + \v aj \ 2 \v pi \ 2 



, n = 3. 



This has an interesting appearance, for it expresses the real part of the box Re 
magnitudes of the four complex V^s which define the box. 



(67) 



in terms of the 



The unitarity constraints ([59]) and (60) greatly simplify our expressions for a doubly-degenerate box m O ai = \ V a i\ 2 : 

J2 v ^a ^rjx) ("E^a ^Ay 



I 1 r 



1 — 1 ( 



^ax oty 

ax n 
<->ay 



i/f/^i, and 



(68) 



(69) 



( Ez#i 

where the first equalities are due to equations ( p^ ) and (p0|). Applying equation (^) to three generations, one finds 
that doubly-degenerate boxes are expressible in terms of the real parts of six ordered boxes, rather than the nine 
complex boxes used in equations (|34|) and ([35]). For example, 



l^lll 4 = =- 



*23 + K33 



71 = 3. 



(70) 
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When considering n > 3, each sum has more terms, but all terms in the numerators in equations fl6q ) and ( p9[ ) always 
contain Rs to the second order, while the denominator terms contain only the first order of Rs. Thus these expressions 
will be much more manageable than equations (|34|) and ( |35| ) which exhibit the fifth order of boxes in the numerator 
and the fourth order in the denominator. 

Summing equation (p9|) over j ^ i yields another expression for |Vju| 2 in terms of nondegenerate boxes: 



i^i 2 E i^i 2 = i^i 2 (i - i^i 2 ) = -EE ( 71 ) 

The same result is obtained by summing equation (^o|) over all j. The explicit solution or the above equation, valid 
for any number of generations, has a two-fold ambiguity: 

1 



V j¥» ^ 



(72) 



When the double sum is small, an approximation to the exact equation (f72|) is 



(73) 



For three generations, this approximation yields \ V\i\ 2 as a linear equation of four nondegenerate boxes. For example, 



\V U 



1 ± sf\ + 4 ( ft 2 2 + &32 + $23 + * 



33 , 



[- ( ^22 + ^32 + ^23 + ^33) , 1 + ( ^22 + ^32 + ^23 + , U = 3. 



C. Unitarity Relationships Relating Real Parts of Boxes with Imaginary Parts of Boxes 



We may use the homogeneous unitarity conditions ( |59| ) arid (|60j) , along with the tautologies (|29|) and ( |30| ) to obtain 
constraints between nondegenerate boxes, thereby reducing the number of real degrees of freedom. Recognizing that 
the tautology (|9|) gives 



™U aj = Re ' "> AJ 



^Aj KAj + ^Aj «Aj - Tr,j KAj JAj + i% ^Aj JAj 



Cra,) 2 + (^) 2 



the unitarity constraint fl59[ ) becomes 
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(74) 



« 'Ha, 'Ha. + Ha- - %i Ha, ^ + ^ J Xj % = - ( ( H . )'' + ( Xv,) ') E ( 75 ) 



with the usual inequalities r\ 7^ A 7^ a, i =/= j satisfied. A similar treatment may be applied to equation (|3C|) to arrive 
at the relation 

°Apx °Ap y "&0y + °?0x °$f3y °&0y ~ °$/3x < R/3j/ "fpy + °&f3x °Sfiy "Sfiy = — (( "R-Py) 2 + ( ^Py) 2 ) °&l3z, (76) 

with a 7^ /3, and x ^= y =/= i. These constraints interrelate imaginary and real parts of boxes through unitarity for any 
number of generations. 

We may infer constraints similar to equations ( f75| ) and (|7^) by setting the imaginary part of a singly-degenerate 
box to zero: 



Im(-D a3 )=Im( XJ )= (77) 



J Aj 

(*4Qrt?M m m m qii m ad ad rid 

•Urjj K\j — "3 n j TUj KAj — T%; T*Aj JAj ~ Tqj gAj JAj _ () 

K 2 Aj + ^ 



(^^q,! 7^ j), and 



cm ori cm; art cti (M m ai a* cwi cwi a* 

^/3x- *\3y ~ ^fix ^py ~ ^(3x^f3y Jf3y ~ Jf3xJf3y ?f3y _q 

(a 7^ /?, and x ^= y ^ i), so 

°^rjj ^Aj ^?7J ^Aj ^Aj ^Aj ^Aj ^r/J °^Xj ^Aj — Oi (^^) 

(r) ^ A 7^ a, i ^ j), and 

oj (8 og is a) a ffl} *} ai oti ad ax n fon\ 

n-px J0y tt/3y — Jf3x K-0y ri0y ~ ^0x *t0 y Jpy ~ J fix Jf3y <>{3y = U. [OV) 

(a 7^ j3, and x ^ y ^ i). With a bit of index switching, these constraints lead to the isolation of a single R: 

OA ad fH \ ad ad 3d, ad ad aA \ ad ad ad 

^ = ggj + ill = ^A? .gft/ + "fl? ^Pv r 8 i) 

with r]^X^= a, in the first expression, and a/ft and x ^ y ^ i in the second. 

D. Indirect Measurement of CP- Violation 

Suppose CP is conserved. Then = for all index choices. The implication for equations ([75]) and (|7^) are 

%i 'Kaj + tix, ^ =0, ( v ^ A + a, and j), (82) 
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and 

%x % y + °Rft, %* = °> (x^y^z, and 0). (83) 

If cither of these relations is violated, then so is CP. An attempt to determine the validity of these relations may be 
facilitated by substituting — ^ 2«=i l^ai| 2 l^aj| 2 f° r the sum over r and — \ | Vq.^ | 2 | Vai| 2 for the sum over z, as 

given in equations ( |59| ) and (|60|). 

For three generations, t must equal A in equation J75|), = — Ta-,- = = »7 by equation (Q, and equation 
may be solved for J 2 directly: 

J 2 = %i % Xj + %i H Xj + <K Aj $ Aj ;, n = 3, (84) 



which expresses ^7 in terms of 3 i?s summed in pairs. Similarly, in equation (|76[), z = y, and = — = by 
equation (|55|), leading to another expression for J' 2 : 

J' 2 = £ K J a J - °&p y + < K/3j + °K/3i/ < R/3y J n = 3. (85) 

These relationships between the Rs and Js in three generations exhibit a simple parameter symmetry: J 2 in equa- 
tion ( |8~I| ) equals "Jt^ terms cyclic in (a, (3, A); in equation (|S5|), it equals 1fyy lft ( 8j-+ terms cyclic in (i,j,y). 
Equations j8~i| ) and ([Sq ) say that the three real elements in any row or column of the matrix B may be summed in 
their three pairwise products to yield the CP-violating invariant J 2 . Note that if CP is conserved and J is zero, then 



equations (84) and ( pig ) also tell us that all three Rs in any column or row cannot have the same sign. 

If CP is violated, J ^ 0, so the combination of real parts on the right-hand sides of equations fl84| ) and Jig 1 ), 
measurable with CP-conserving averaged neutrino oscillations, cannot be zero. Thus, even if CP violating asymmetries 
are not directly observable in an experiment, the effects of CP violation may be seen through the relationships 
among the real parts of different boxes! A related observation has been made for the quark sector by Hamzaoui in 
reference jL9|. 

For three generations, unitarity enforces the simple relations (jHJ) and (|5^) among the Js in equation ([u]), which 



then reproduces the expressions ( p4|) and (|85|). For n > 3, however, equation (plf) is different from the latter equations 



since the sums ^2 T ^ a b anc ^ Xz^ij m equations (75) and ( |7q ) ensure that boxes with all n flavor or mass indices 



enter into the equations replacing (p4[) or (|85j), respectively. Equation (81) involves only three such indices for any n, 
so it will be different from the unitarity constraints in higher generations, providing additional identities among the 
real parts of boxes. 
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E. Inhomogeneous Unitarity Constraints 



The homogeneous unitarity constraints cannot provide the desired normalization of the V a i or the boxes; we need 
the inhomogeneous unitarity constraints arising when the Kronecker delta in equations ( ffO| ) and ( ffl| ) is not zero. We 
now develop these inhomogeneous constraints^ which are functions of degenerate boxes and therefore purely real: 



m U m + ma m = V mn «* - and 



(87) 



(XI 



E 



a/ 



-/ cein 



Parenthetically, we note by comparing equations (|87]) and ( |88| ) that a sum over mass-degenerate boxes equals a sum 
over flavor-degenerate boxes: 



(89) 



Equations (|87j) and (|88|) can be rewritten strictly in terms of nondegenerate boxes by using the homogeneous 
unitarity constraints ( |59| ) and ([30|). We hnd 



E -E%*]= ' 



(90) 



with x ^= y ^ i, and 

J2 x ^i °^Ai) (- E 



y=£i n vi 



( Et/i ftrii 



\ 



(-E,^ ^a.) (-E 



E -E°M = ' 



(91) 



Et^i ^ni 

with \ =/= rj =/= a. Note that these inhomogeneous unitarity constraints do not involve the Js. 

Isolating the square root, squaring the equation, and multiplying through by the resulting denominator, we get 
quartic e quations, each relating n(n — 1) Rs: 

3 In terms of mixing-matrix elements, equations (Jst|) and ( |S8| ) are trivial. For example, equation (^) is 

n 

| , M 4 +^|V^| 2 |1^| 2 = | , M 2 , or ^2\Vm\ 2 = 1. (86) 



'1=1 
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1 + 2 

2 



(92) 



(x ^ y i), and 

for A 7^ 77 7^ a. In three generations, each of the sums in equations (|9^) and (^3|) has only two terms, which is not so 
formidable. For example, for a — 2 and i = 1, we have 



(93) 



;£ 22 + % 32 )(£ 23 + : ft33)(% 



23 



V33 



[1 + 2 ( % 22 + 1 3 2 + %3 + % 3 )] + 



and 



( %2 + % 32 ) 2 ( %3 + %33) 2 + ( &22 + %2 + % 23 + %33) 2 ( $23 + ^33) 2 = 0, U = 3, 

' ^ 22 + ^ 23 ) ( %2 + ^33) ( ^32 + & 33 ) [1 + 2 ( &22 + %3 + %2 + %3 3 )] + 
( % 2 + ^ 23 ) 2 ( % 32 + ^SS) 2 + ( &22 + ^ 23 + %3 2 + ^SS) 2 ( %2 + &3 3 ) 2 = °' 71 = 3 - 



(94) 



(95) 



V. REDUCTION TO A BASIS 



A. General Algorithm 



As with the constraints on the imaginary parts of boxes, many of the constraints ( |74| ) to (93) are redundant, but an 
independent set may be used to reduce the number of box parameters to a basis. The homogeneous constraints (|75|), 
(fZ6|) , and ( |8l| ) are much simpler than the inhomogeneous constraints (|9|) and (|93|), so it is advantageous to use as 
many homogeneous constraints as possible to construct the basis. Still, some inhomogeneous constraints must be 
invoked if the boxes and matrix elements of V are to be normalized. 

Notice that with only fixed i and j mass indices and a flavor sum over the entire column of B, in equation (|75|), 
and only fixed a and (3 flavor indices and a mass sum over the entire (a, (3) row of B in equation (|76|), these equations 
express relationships among all the boxes within single columns and rows, respectively, of the matrix B defined in 
equation (EJ). For each choice for the set of indices, equations (ra) and (|7q) may be used to express one R or J in 



terms of the others. This suggests a general algorithm for constructing a basis set of boxes. One R from each row 



and column may be taken as dependent using the homogeneous unitarity constraints ( |75| ) and (|76|), or (81). This 
leaves (n — l) 2 Rs and q Js in the independent set, where q is the number of Js which are still independent after 
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the application of the homogeneous constraints equations ( pij ) and ( p5| ) . The n x n mixing matrix is parameterized 
by (n — l) 2 parameters, so we seek a basis set of (n — l) 2 box parameters. The inhomogeneous unitarity equations, 
which we will address shortly, must supply the remaining constraints. The number of inhomogeneous constraints to 
be used is apparently equal to q. We have seen that with three flavors this is just one. The end result is a basis of 
(n — l) 2 Rs and Js with no more than n — 1 Rs coming from any row or any column. If one restricts the bases and 
counting to accessible boxes, the algorithm may differ for n > 3. We do not pursue this here. 



B. Three Generations 



We provide here an example of a basis construction for three generations obtained by substituting in the unitarity 
equations derived above. Rearranging the three-generation equations (|84| ) and ( |S5l) yields expressions for one R in 
terms of two other Rs and J: 

We may then eliminate ¥£23 by either equation (Bq) or equation (B7|) 



12 _ ^33 ^33 ~ J 2 _ ^22 ^23 - J 2 / qs x 

23 ~~ - ^33 - K33 ~ - S.23 - &22 ' [ j 



We may similarly eliminate ^32 and ^,33 



u ^22 ^32 — J 2 _ ^33 ^33 ~ J 2 /qqn 



^32 - ^22 ~ $33 - $33 

and 

^23 $33 — J 2 $32 $33 — J 2 



24 ^23 ^33 — U ^32 J^33 — U ,. nn , 

%3 ~ ~ $33 - %23 " -^33-^32- (100) 

$33 may be eliminated from equation (pq), and $33 may be eliminated from equation (|100|): 



u $33 $22 $23 + ( $23 + $22 " $ 33 ) , m , 

$33 $23 +^$22- $22 $23+^ ' ^ ^ 

U _ $33 $23 $32 + J 2 ( $32 + $33 ~ $2 3 ) . . 

33 $23 $32+ ^33 $23- $32 3*33 + J 2 ' 1 j 
Equation ( |99| ) will not provide an additional constraint; it is redundant to the other two. As expected, we must turn 

to the inhomogeneous constraints to eliminate the last degree of freedom. We will here choose the constraint (|94|), 
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since its expression in terms of our four remaining boxes is the least complicated. We may substitute the second 



t 3 s)+J z 

|2 



equality from equation (IOC) for ^33 and the second equality of equation (|9^) for ^23 into equation (|94|). 

Multiplying through to place all of the terms in the numerator, we are left with a constraint which contains the 
five parameters $22, %32, ^33, and J 2 : 

= ( %2 + % 3 f ( ^22 + 1 32 ) 2 [- %2 % 3 + %3 ( %2 + I33) + J 2 } 2 
%32) 2 + ( ^23 + + ^22) (132 

X [ % 3 ^33 + ^22 (- ^23 + ^3 3 ) J 2 } * (103) 
^22 + 1 32 ) ( %2 + % 3 ) [ %2 ( ^23 ~ 133) - ^23 1.33 ^ J'] 

%2 + 2 £22 I32 + 2 ( 1 32 ) 2 + 2 % 2 % 3 

+ % 3 + 2 % 2 133 + 2 % 3 % 3 + 2 J 2 ] [ % 3 ( % 2 + I33) ~ %2 133 + J 2 ] ■ 

This equation is quartic in all five parameters. We may eliminate any one by either algebraic or numeric means, 
leaving us with the desired four parameters as the basis. 



VI. SUMMARY 

Neutrino physics has entered a golden age of research. New experiments all over the globe promise an unequaled 
amount of data from the sun, the atmosphere, accelerators, supernovae, and other cosmic sources. The latest data 
suggests that more than three neutrino flavors may participate in neutrino oscillations Q. Analyzing such refined 
data requires a consistent, model-independent approach which may be easily applied, and easily extended to higher 
generations. Our work presented here offers such an approach. One result which we view as particularly noteworthy 
is that high-statistics data on averaged oscillations is sufficient to determine the conservation or non-conservation of 
CP in the lepton mixing matrix; this indirect test of CP is a consequence of the unitarity of the mixing matrix. 
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APPENDIX A: GRAPHICAL REPRESENTATION OF BOX RELATIONS 



Many of the relationships between boxes developed in Section III E can be derived using a graphical representation. 
In this method, boxes in the numerator of a product are represented by two vertical lines; boxes in the denominator are 
represented by two horizontal lines. Lines exit the locations of the matrix elements which are not complex-conjugated 
and enter the locations of the complex-conjugated matrix elements. For example, the box 11 D 22 = VnVft^V^V^* 
is represented by a vertical line pointing from Vu to V21 and a vertical line pointing from V22 to V12, as shown in 
Figure |a. The inverse box 11^ = ( 11 ^22 ) is represented by a horizontal line pointing from Vu to V12 and one 
pointing from V22 to V21, as shown in Figure ^|b. The complex-conjugated box 11 UT 22 is equal to 12 U 2 \ , so one just 
reverses the arrows to complex conjugate a box, as shown in Figure |^c. 



(a) (b) (c) 

FIG. 5. The graphical representation for (a) n Q 2 2 , (b) ( n Q22 )~ , and (c) u n| 2 . 

To multiply boxes together graphically, one merely draws the lines corresponding to each factor on the same grid. 
Horizontal arrows entering or leaving a point cancel out vertical arrows entering or leaving, respectively, that point. 
Uncanceled arrows entering a point signify the survival of the complex-conjugated matrix element associated with that 
point. Those leaving a point signify the survival of the ordinary matrix element of that point. Figure ^a represents 
the product n D22 12 ^3Z ( 23 32 ) ■ The horizontal arrows do not cancel vertical arrows at any point, so no 
simplification may occur. The upward arrow at V22 represents that element in the numerator. The horizontal arrow 
represents the element V 2 \ in the denominator, which does not cancel. Counting off the arrows at each vertex, we 
find the expression 

"□ 22 12 n 33 ( 23 n 32 T 1 = VllVl2 ^f 3 ^f" 2V3 \ (Al) 

V 22 V23 V 32 V33 

which agrees with the definitions of boxes. 
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(a) (b) 
FIG. 6. The graphical representation for the products a) n a 2 2 12 a 33 ( 23 D 32 ) ~\ and b) n D 2 2 12 Q33 ( 22 Q33) _1 - 

Figure ||b represents 11 22 12 ^33 ( 22 ^33 ) , a product in which some canceling does occur. Picking out the 
uncanceled arrows at each vertex, we are left with 



11 12 / 2 2 n \-l _ ^11^12^1*2^*3^* 
u 22 u 33 ^ u 33 ) — 7JZ 

The graphical method is a powerful tool for finding relationships among boxes. For example, trying to obtain the 



□22 ( "Daa ) _1 = (A2) 



equations (34) and (Ba) for 



m O ai = iVai] 4 of Section HIE without graphs involved quite a few false starts. Using the graphical method, we 
need only find a series of arrows which cancel for every point except (a, i) and leave two incoming and two outgoing 
vertical arrows at that point. Consider | V21 1 4 as an example. We choose to draw all of the arrows involving V21 
pointing downward, as shown in Figure 0a. These arrows must be part of boxes, so in Figure we add the arrows 
to finish those boxes. Next we draw two horizontal boxes in Figure ^c to cancel the extra arrows in the first column 
of the matrix. This still leaves V22 and V23 with two sets of uncanceled arrows apiece. In Figure [?]d we draw two 
more horizontal boxes to compensate. This adds arrows to our previously clean V12, V13, V32, and V33. Drawing the 
final vertical box in Figure 0e cancels those. Recapping what we have done, we see that step (b) completes 11 CI22 , 
11 D 2 3 , 21 CI32 , and 21 CI33 in the numerator. Step (c) divides by 11 CI32 and 11 CI33 , and step (d) divides by 12 D 2 3 
and 22 ^33 . Step (e) multiplies by 12 CI33 , leaving only the point V21 with uncanceled arrows. It has two vertical 
arrows coming in and two leaving, so our graph represents the equation 

21 n 11 n 21 n Hn 12 n 

i T r 14 u 32 u 22 u 33 u 23 u 33 

\V21\ - 



11D 32 11D 3 3 22 D33 12 D23 



of Section III E . Other examples of this representation are included in that section. 
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(d) 



(e) 



FIG. 7. The steps to obtaining | V21 1 4 as a function of ordered, non-degenerate boxes. The additions in each step are 
designated by the thicker arrows. 
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